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Abstract 

We calculate the beta function of non-linear sigma models with S D+1 and AdSo+i target 
spaces in a 1/D expansion up to order 1/D 2 and to all orders in a'. This beta function 
encodes partial information about the spacetime effective action for the heterotic string to 
all orders in a'. We argue that a zero of the beta function, corresponding to a worldsheet 
CFT with AdSo+i target space, arises from competition between the one-loop and higher- 
loop terms, similarly to the bosonic and supersymmetric cases studied previously in pp. 
Various critical exponents of the non-linear sigma model are calculated, and checks of the 
calculation are presented. 
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1 Introduction 

Particular interest attaches to backgrounds of string theory involving AdS^+i because of 
their relation to conformal field theories in D dimensions [3E1I1] (for a review see 0). But 
because these geometries (with some exceptions) arise from the near-horizon geometry of 
D-branes, formulating a closed string description is complicated by the presence of Ramond- 
Ramond fields. 

It was recently proposed £Q that AdSp+i vacua might exist without any matter fields at 
all. Instead of relying upon the stress-energy of matter fields to curve space, the proposal 
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is that higher powers of the curvature compete with the Einstein-Hilbert term to produce 
string-scale AdS^+i backgrounds. The main support for this proposal comes from large D 
computations of the beta function for the quantum field theory on the string worldsheet. 
Before discussing these computations, let us review the lowest-order corrections to the beta 
function in an a' expansion: 

bosonic: 
heterotic: 
type II: /% 

These expressions are obtained using dimensional regularization with minimal subtraction, 
and all derivatives of curvature are assumed to vanish as well as all matter fields. Derivatives 
of curvature indeed vanish for symmetric spaces: for example, 

Rijkl = —J2^9ik9jl — 9il9jk) (2) 

in the case of AdSn+x- One indeed finds non-trivial zeroes for AdSp+i from all three beta 
functions in (JIJ). An examination of higher order corrections in the bosonic and type II cases 
shows that the zero persists in the most accurate expressions for the beta function that are 
available at present; however its location changes significantly, converging to a'D/L 2 = 1 as 
D becomes large. One aim of the present paper is to pursue similar large D computations 
in the heterotic case. 

It should be clear from the outset that the question of the existence of AdSn+i vacua with 
a'D/L 2 close to unity is a difficult one to settle perturbatively. Fixed order computations 
are not reliable guides because the scale of curvature is close to the string scale. Large D 
computations with finite a' D / L 2 seem to be a better guide, but they too could be misleading, 
mainly because higher order effects in 1/D than we are able to compute could change the 
behavior of the beta function significantly. These difficulties were discussed at some length in 
PP. Also, the vanishing of a beta function such as the ones in is only a necessary condition 
for constructing a string theory: one must also cancel the Weyl anomaly and formulate a 
GSO projection that ensures modular invariance and the stability of the vacuum. 

There is a more general reason to be interested in high-order computations of the beta 
function on symmetric spaces: from them we can extract information about the structure of 
high powers of the curvature that is quite different from what is available from expansions of 
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= a'R l3 + —R tklm R 3 klm + 0(a' 3 ) 
a' 2 

= a'Rij + —R lklm R 3 klm + 0(a' 3 ) (1) 

= a'Rij + — R m hkiRjrt m {R k qs r R tqSH + ' qs* R hrS9 ) + 0(o/ 5 ) . 



the Virasoro amplitude. While the latter tells us about terms involving many derivatives but 
only four powers of the curvature (because only four gravitons are involved in the collision), 
the former tells us about many powers of the curvature with no extra derivatives. 

The organization of the rest of this paper is as follows. In section |21 some general 
properties of the heterotic NLaM are discussed. In section |3J the formalism and the results 
at 1/D order are presented. In section the critical exponents at 1/D 2 , the beta function, 
and the central charge of the CFT are computed. The appendices include a brief explanation 
of the method of the calculation for the diagrams needed and the values of these diagrams. 



2 The heterotic non-linear sigma model 

As in P , much will be made of a connection through analytic continuation of the NLcrM on 
AdSo+i and the NLcrM on S D+1 . If L is the radius of S D+1 and g = a' /L 2 , then continuing 
to negative g leads to the AdSo+i NLcrM. The argument in |]Q is slightly formal because it 
relies on an order-by-order perturbative evaluation of the partition function. 
The action for the S D+l heterotic NLcrM is 

S = / d 2 xd9 LD + $0_$ + A($ 2 - 1)1 + — f d 2 x\ A d + \ A (3) 
Aug J Ang J 

where 

$ = S + 9* A = u + 9(r D+ = -= + 9— d ± = —. (4) 

t)9 OX-\- OXzf: 

A is a spinorial superfield, and u and \& have opposite chirality. This leads to the action 

S = — [ d 2 x \(dS) 2 + Md^ + crfS 2 - 1) + 2u<5>S] . (5) 
4irg J 

We have omitted the fermions from (jSJ) because they decouple from the gravitational 
action when the gauge field is set to zero [01 E| as in our case. The Feynman rules for the 
theory (jJJJ) can be seen in figure There is also a tadpole for a, but we omit it because it 
does not contribute to the Dyson equations for the scaling parts of the dressed propagators, 
as in [S]. 

After a change of variables that renders the kinetic terms canonical, we can continue 
to negative values of g as in pQ to obtain an AdSo+i heterotic NLcrM. Quantities that are 
computed locally and perturbatively, such as n-point functions, cannot distinguish between a 
space of positive or negative curvature. As the beta function is derived from such quantities, 
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Figure 1: The Feynman rules for the heterotic sigma model. The shaded circles indicate a 
dressed propagator. The circles indicate that a loop involving only the components of the $ 
superfield receives a factor of N. We have suppressed the tensor structure of the rules since 
only 5^ appears. 



it too can be continued to negative g, at least order by order in perturbation theory. 

The heterotic NLaM on S D+1 is a generalization of the 0(D + 2) model, and much of 
the relevant literature concentrates on an expansion in 1/(D + 2) rather than 1/D. We will 
therefore set 

N = D + 2 (6) 
and work with N or D, according to convenience, in the rest of this paper. 



2.1 Some properties of the heterotic NLcrM for large D 

It is known that in the bosonic sigma model a mass appears jH] in the 1/N expansion. The 
same phenomenon appears in the supersymmetric extension of the sigma model fOJ ITT] 
where also the fermions acquire the same mass, signaling chiral symmetry breaking. In the 
heterotic case the bosons S also acquire the same mass, showing that the interaction term 
does not destroy this effect. To understand this, let's start from our action (j3J, and in the 
partition function integrate first the fermionic fields and then the bosons, since the action is 
quadratic in these. We have omitted normalization factors of the partition function in the 



following. 



VSVmVaVu exp 



Aug 



J VSVaVu [det^J^exp 

/ 



d 2 x[{dS) 2 + + a(S 2 - 1) + 2u(S ■ *)] 



4- [ d 2 x (s(d 2 - a)S - S^uS'' 



VaVu [det(z^)] 
/ VaVu e iSaS , 



det ( — d 2 — au—u 



exp 



Aug 



d xa 



(7) 



where the effective action for the Lagrange multiplier fields is given by 
S P fr = I d 2 x 



Y—cr - — Tr log(z^) + — Tr log ( -d 2 - a - u^-u 
Aug A 2 \ lip 



Since we are taking the limit N — > oo with g N finite, we see that all terms in the action are 
of order N. We can evaluate this integral by the method of steepest descent, i.e. by finding 
the classical value of a, u that minimizes the exponent, as is done for instance in [T21 113j . 
This gives the variational equations 
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{x 



-d 2 — a(x) — u-hU 



=-j— |z) 



1 



2nNg 



(9) 



-d 2 — <j(x) — u^rKi 

V 



=T-|a;) = 0. 



Because the right hand sides are constant, the left hand sides must also be constant. A 
solution to these equations is given by 



u{x) = cr(x) 



-m 



(10) 



It is easy to see that j^u(x) = const, has as its only solution u = 0. This is in contrast to 
the super symmetric case [T2], where there are three solutions. Now m 2 must satisfy 



d 2 k 



1 



1 



k 2 + m? 2i\g§ 
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Using a simple-momentum cutoff, 2v ] joN = 2^ m • By renormalizing at a scale M we get 

1 1 1 , A 

^ lo ST7- ( 12 ) 



2%g N 2ngN 2n M 

Solving for the mass m we get m = M exp[—l / gN]. Since this is a physical mass we expect 
that it does not depend on the renormalization scale. Using the Callan-Symanzik equation 
for m, 

M ^T7+^W }m(g,M)=0, (13) 



DM ^^ J dg 

gives the beta function (3(g) = —g 2 N. The mass is the same in the bosonic, supersymmetric, 
E] and heterotic model. This could have been predicted since the first order (3 function 
is the same in all models, (3^ = a'Rij. Another way to get the same result is to calculate 
the effective potential for the a field and see that the minimum of the potential is not at 
zero but at a = Me~ 1/9N . One can go further and examine the effective action (JHJ). It is 
easy to evaluate the counterterms needed for one-loop renormalization, as we have already 
computed the wave function renormalization of the a field, and doing so one finds 

£ , eff = J-^l + gN\og^ja-j Tr log + | Tr log (-9 2 - a - fi^u) . (14) 

The bare and the dressed quantities are related by 

a = Za u = Z 1 / 2 u 9o = Z- 1 g Z = 1 + ^ log ^ . (15) 

Calculating the quadratic terms in the fields u, a will give us the propagators for these fields. 
We can easily find 

q 1 _ 1 N 1 1 

S e s = — Tr— - -u— u- — Tr— 2 a ^2 2 a • ^ lb > 

2 —o z — m z zip — m A — o z — m z —o z — m z 

The next step is to evaluate the propagators. One finds fTJ 

2? 2? 
S u (k) = --(# - 2m)V(k 2 ) D a (k) = -(Am 2 - k 2 )V(k 2 ) , (17) 

where in d dimensions 

(AttY/ 2 /Am 2 - k 2 \ 1 ~ d/2 ( k 2 
V(k 2 ) = ( A ^—l- L Fl (2 - d/2, 1/2, 3/2; 



AT(2-d/2)\ 4 J V ' ' 'k 2 -Am 2 
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In two dimensions this simplifies to 

^='E(- iih 65) • (i9) 

It is easy to see that the u propagator in d = 2 dimensions has a pole at k 2 = 4m 2 . This 
means that there is a particle with this mass. Since the classical equations give u = —iS<ft^, 
there is a boson-fermion bound state, created by the operator S^. This is in agreement 
with [H], but there is no supersymmetric corresponding fermion-fermion bound state, since 
the Gross- Neveu interaction that is responsible for it is absent. 

One can also extend the calculation of [T2*l HI] to show that there is no multi-particle 
production in the heterotic sigma model. For example, the process 2 — > 4 particles can be 
shown to vanish. The reasoning is that the formalism of [T3|, valid for the bosonic case, can 
be extended to include superfields. 



3 Critical exponents in the 1/D expansion 
3.1 General discussion 

The method used to determine the critical exponents is the one developed in jHJ for the 
bosonic model and extended to the supersymmetric case in ^H] Ej- To this end one writes 
expressions for the propagators of the fields near the critical point. In keeping with the 
notation of (HI H3 EE] we assign dimensions to the fields 

dim[5] = dim[tf] -- = A s = (d-2 + rf)/2 

I (20) 
dim[a] = dimfw] + - = A CT = 2 — rj — x ■ 

For small but non-zero x, the two-point functions may be expanded as follows: 
°ss(x) = -^-(1 + TgsX ) GW(z)= ^ (l + r w ^ ) 

X Z X ^ 

G„ a {x) = %(1 + K a x 2 ") G uu (x) = ^^^(1 + 2A ) , 

where 

7p = PV _M (22) 
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is the chirality matrix in 2 dimensions. We have omitted the O(N) indices because both 
the propagators and the vertex are proportional to 5^. So all Green's functions can be 
expressed as a scalar function times products of 5^, and one does not have to worry about 
tensor structures of the form (x — y)^{x — y) u . Then one can write the Dyson equations in 
a 1/D expansion for the propagators. Graphical expressions of these equations are shown 
in figures 121 through |3| The Dyson equations impose consistency conditions on the critical 
exponents that determine them completely. The graphs that appear in the Dyson equations 
are the 1PI graphs with the exception of graphs which contain subgraphs that already 
appear in the Dyson equations at a lower effective loop order: in other words, we exclude 
diagrams that are already taken into account by expressions for the corrected propagators. 
The effective loop order is the number of loops minus the number of loops involving only the 
components of the $ superfield. 

The left hand side of each Dyson equation is a 1PI propagator, which is the inverse of 
the connected two-point function. These inverse propagators are computed by first passing 
to Fourier space using 

V = A, • < 23 > 

The inverse propagators are found to be 1 

G~Hx)= l+lP r ^-Vt (l-s(A -1 X)T' x 2X ) 

* 1 uu-L 



(24) 



where 



H = d/2 A s = fi-l + r]/2 A a = 2-7]-x, (25) 



J Note that Gq,^ ■ G^,\, does not strictly give the unit matrix but instead 1 ~t 2 7p = f J j in two 

dimensions, which is what we want. Also note that in the Dyson equation for \I/ in the right hand side one 
encounters u propagators that give the right chiral structure, and vice versa for the Dyson equation of the 
u field. 



s 



and for arbitrary y, 

a(y) = r(/i = V) p(y) = a{V ~^ r(y) = ^ 

a{y- X)a(y + X-fj.) , . y(y - fj,)q(y, A) 
« fe ' A) = afeW,-,) «(»■*)= (1/ _A) (1/ + A 

To calculate the beta function, one first evaluates the critical exponents of the model at the 
fixed point in d = 2 + e dimensions and then uses the relation 

A = - l -(3'{g c ) , (27) 

valid at the critical point, to extract (3(g). This is possible because, in dimensional regular- 
ization with minimal subtraction, the only e dependence in (3(g) is an overall additive term: 
see £Q for details. Expanding in 1/D with k = gD held fixed, one finds 



A(£) = E MfJ /M. £ _ K + E w (28) 



Ai(e) (3(g) _ , 

i=0 y i=l 

where 



Note that the critical exponent A is a measurable quantity and as such it should not depend 
on the renormalization scheme used. Passing from the Aj(«) to the bi(n) does introduce 
significant scheme dependence. 

3.2 Critical exponents at order 1/D 

The Dyson equations can be expressed in terms of parameters 

which can be regarded as dressed vertex factors for the two vertices shown in figure H The 
leading non-trivial Dyson equations come from the graphs labeled So, a, S 0) b, $o> n , and 
F in figures 121 through |3J the tree-level graphs make no contribution to the leading scaling 
behavior. The quantities T' ss , T'^, T' a(T , and Y' uu describe how far one is removed from the 
fixed point; so in particular one must be able to set them all to zero and get a self-consistent 
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Figure 2: The Dyson equations for the S propagator. 

set of equations. Then the dependence of each graph on the position-space separation x is 
just an overall power of x. Matching these overall powers leads simply to the constraint 
X = 0. Matching other factors leads to the equations 



p(A s )+w + v = r(A s )+v = 

^r(A a - 1) + v = ^p(A a )+w = 



(31) 



which determine the quantities As, A a , w, and v as functions of \x and N. The system is 
in fact over-determined if we recall the relations ()25|) and the constraint x — 0. But we will 
see in section E~T1 that x = is only a leading order result; thus to solve (|31|) we expand 



i>0 



— W = > : V = > : 

i>0 i>0 i>0 



(32) 



Then one straightforwardly extracts from (|31|) the coefficients 



^70 = xo = 

m = -2 



w = 

r(2//-i) 



v = 



(/i-ifrii-^-i^+i; 
(2-/x)r(^-i)r(/i + i) 

Wl = 2^ * 

;i-/i)r(/i-i)r(/i + i; 



2tt 2 ^ 



m ■ 



(33) 
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Figure 3: The Dyson equations for the ^ propagator. 




Figure 4: The Dyson equations for the a propator. 

Higher order coefficients receive contributions from higher order graphs. Note that Xo = 
could be obtained either from matching overall powers of x or from the equations (|31j). 

Now consider non-zero coefficients T' ss , T^, T' acr , and T' uu : this corresponds to moving 
away from the fixed point. Linearizing the Dyson equations leads to the constraints 

-p(A s )q(A s , x)V ss + w(T' ss + r' aa ) + «(ru + r;j 

-r(A s )s(A s ,X)r^ + v(T' ss + r uu ) 
-r(A CT - l)s(A a - 1, X)T' UU + Nv(T' ss + r' w ) 

- P (A CT )g(A CT ,A)r; (T + ivr' s ^ 

Graphically, these equations arise from using the leading power-law expressions (e.g. Tss/x' 






0. 



(34) 
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Figure 5: The Dyson equations for the u propagator. 

rather than Gss{x)) for all propagators except one, chosen arbitrarily; and for that one, use 
the correction term (e.g. r ss^'ssl 'x 2As ~ 2X )- The linear equations (J3*4^ must admit a non-zero 
solution for T' ss , T^, T' aa , T' uu in order for the correction terms to describe a genuine de- 
formation of the critical point. So the corresponding determinant must vanish, which leads 
to 



O + v)q(A s , A) + w 1 



q(A a ,X)J s(A a -l,X) 
l-s(A CT -l,A)) 2 



1- S (A CT ,A) S (A CT -1,A)) 

(35) 



s(A CT -l,A) 



Note that setting v = gives the equation valid for the bosonic model as expected. To 
simplify ()33|) . one can use and note that W\ = V\j^. So far we have not used any 
expansion in 1/D. Using the expansions ()32|) and 

A = E4 (36) 



we can determine 



D i 

i>0 



A = //-l A 1 = i(2/i-l)(^- 1)771. (37) 



Another way to compute the original determinant, is the following: one notes that r(As)s(As, A) 
1/D°, while all other terms scale at least as 1/D 1 . This means that, in expanding the 4x4 
determinant of the first order contribution comes only from the determinant of the 
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3x3 matrix 



AI 



V 



-p(A s )q(A s ,X)+w 
w 

V 



-p(A a )q(A a ,X)/N 




v 




(3J 



-r(A CT -l)s(A (J -l ) A)/iV J 



where in each element of the matrix we only keep the first term in the 1/D expansion. This 
determinant provides the first 1/D term of (|35j) and thus reproduces (|37[) . 

It is interesting to compare with the bosonic case. It is easily seen that A^ 10 * = |A^ os , 
even though rft et , wi, and V\ are not so simply related to ri\ os and the corresponding vertex 
factor for the bosonic case. The relation X\ ct = ^X\ os is expected: we know that A™ p = in 
the type II case, and having half the fermions in the heterotic case will cancel only half the 
bosonic contribution. 



3.3 A check of the calculation 

Following jTHj, we see that rj is the anomalous dimension of the S propagator. So far we 
have computed it in the 1/D expansion using techniques in position space. One can also 
straightforwardly compute rji in momentum space using the expressions for the propagators 
that we have previously found ()17j) . Firstly we note that for small k, 

Gss(k) ~ £r 2+ " ~ k~ 2+r >° (l + ^ log k + 0(1/ N 2 )) , (39) 

with Gss the Fourier transform of the S propagator. But Gss{k) can also be determined 
from the one-particle irreducible diagrams E(/c 2 ): 

Gg&h) = k 2 + E(k 2 ) - S(0) ~ k 2 (l - ^ tog*) , (40) 

since 770 = 0. Having calculated the propagators of the lagrange multiplier fields it is straight- 
forward to compute the S(/c 2 ) from the one- loop diagrams. We find 

nk 2 ) = 1 * A.CP 2 ) - f 7^t r7 -i-^ - - 2 f ddp timv{p 



(2ir) d (p - k) 2 u ^ ' J (2n) d $ - f) uvn NJ (2vr) d (p + k) 2 

d ~ 2 /NJ^ f M , r in p 2 k cos6sm d ~ 2 

dp / dO- 



2 Fi(2-d/2,l/2,3/2,l)r(2-d/2)r(^i) J ^ J p 2 + k 2 + 2pk cos 9 

(41) 
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where we have put m? = as usual [TZj, and M is the cutoff. One notes that, for small k 
the k 2 \ogk behavior comes from the small p region [TBj . The integral is trivial to do, and 
for the k 2 log k part it gives 

E(k 2 ) = -- -rk 2 \ogk. (42) 

1 ; #^(2-^/2, 1/2, 3/2, l)r(2-f)r(f) S 1 ; 

It is easy to see using d = 2/x and properties of the Gamma function that this coincides with 
the expression for 771 in (J33)) . An easier way to do the checking is the following. One can 
write similar expressions for E(fc 2 ) for the bosonic and supersymmetric models [HI El- Then 
it is easy to observe that 

£ bos (A; 2 ) - £ bos (0) = -S sup (A; 2 ) + 2£ hct (£; 2 ) . (43) 



This easily gives 2 
With E7] 



27?het = Vbos + Vsup ■ (44) 



_ (2-/i) 4 r(2/i-2) 

^bos ??sup Vsup -\TT2/ \Tf \ \^) 



IX ,bup " up NT 2 (i2-l)T(2-fi)T(fi) 
we find ?7h e t = -Vsup which agrees with (j3~3"l) . 



4 Results at order 1/D 2 

Each graph in figures |2E carries an overall factor 1/D M where M is the number loops minus 
the number of loops containing only S and To see this, first note that each propagator 
Gxx carries a factor Txx (where X = S, a, or u). Next note that the amplitude for each 
graph must contain an overall factor which is a product of the factors w = r 2 ^ s T rTrT /(27rg) 2 
and v = Tss^^^uu/i^g) 2 , one for each vertex in the graph. The overall factor 1/D M 
arises because w and v scale as 1/D and because each loop containing only S and \l/ carries 
a factor of N. The graphs in figures El and 01 are those with M < 2, and the ones in figures 0] 
and El are those with M < 1. Together, these are all the graphs that can contribute to 77, 
w, v, and A through order 1/D 2 , and they also determine x through order 1/D. Because we 
quote final results in terms of 1/D, we must keep in mind the relation between expansions 
2 We have used that E sup (0) = and S hct (0) = 0. 
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in 1/N and l/D: 



i>0 i>0 



with similar relations for other quantities. 



4.1 Calculation of 7/2 

A technical complication arises in the l/D 2 corrections to the Dyson equations that was 
explained and resolved in [HJ • The problem is that the higher-loop graphs diverge when 
X = 0. In fact, x = only up to l/D corrections. But it convenient to regularize the 
"divergence" and extract finite expressions for the two-loop Dyson equations through the 
following steps: 

1. Shift x ^ x + A. 

2. Expand the amplitudes for individual graphs in powers of A. 

3. Cancel 1/A terms against certain counter-terms in the lagrangian. 

4. Fix xi by setting to zero certain terms in the Dyson equation which depend loga- 
rithmically on the position-space separation x and which, if non-zero, would spoil 
self-consistency. 

We will now go through these steps in detail for the \& propagator. The reader who wishes 
to bypass the technical details can skip to f|56[) and (|57|). which are the two-loop Dyson equa- 
tions with all divergences removed. But the results (|54|) and (|58|) for xi provide important 
consistency checks. 

When x 7^ 0, 1 dependence cannot be canceled out of the Dyson equations in a simple 
way: setting T' ss = T'^ = T' aa = Y' uu = 0, one obtains for the ^ propagator's Dyson 
equation 

r(A s ) + v{x 2 ) x + t» 2 (x 2 ) 2x $i + iYW(x 2 ) 3x $ 2 = . (47) 

Here $1 and $2 are functions of Ag, A CT , and \i which diverge when 2Ag + A CT — 2/j = — x — 0. 
Although these are in some sense an artifact of a limit (% — ► 0) which one cannot take 
independently of the large N limit, it is convenient nevertheless to regulate them, as explained 
above, by shifting 

X ^ X + A. (48) 
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The amplitudes $i 2 may then be expanded as 

= ^ + % + 0(A) , (49) 

where both Xi and are functions of Ag, A CT , and /i, subject to 2Ag + A CT — 2/i = 0. In 
appendix we exhibit $12 in the form (|49|) . as well as a number of related quantities that 
enter into other Dyson equations. To cancel the divergent 1/A terms in $1,2, one may rescale 
the lagrange multiplier fields in the original action (J3J). This rescaling amounts to adding 
counter-terms to the action, and it can be expressed, to the relevant order, as 

v^(l + j£)v w^(l + ^w. (50) 

(The factor on v and w is the same because of supersymmetry.) Subjecting (|4"7)l to the shift 
()48j) and the rescaling (|5L)|). it becomes, keeping terms up to 1/N 2 , 



r(A s ) + (x 2 ) x (v + v 2 <f>[ + Nv 2 w& 2 ) + 



(51) 







The last line contains all the divergent pieces. Setting \ = [Hj and taking the limit A — > 
determines mi as 

m l 2^1 , AT 2 ^2 /r N 

" Wl AT = t,1 A" lWl A"' ^ ' 

Plugging back into (|51[). and now considering a finite \ we get 



r(A<?) + (x 2 ) x [y + ^ 2 $i + Nv 2 w& 2 ) + t^Xx 



^2N2 X _ (x 2 )x 



X 

x 2 ) 3 * - {x 2 Y 



+ Nv{w x X 2 [ x ) = 



(53) 



When one expands \ = Xi/N + 0(N 2 ), there are terms that behave as logx 2 . One gets rid 
of these if xi obeys 

Xi = -viX x - 2v 1 w 1 X 2 . (54) 

We could have derived (|52|) . (|54j) purely within the N — > 00 limit. In this setup there is no 
need for A and x is taken to be finite. Note that it behaves as x ~ 1/-^ since Xo = (|33|1 . 
The Dyson equation after the rescaling (jBTIj) is (JBTj) with A replaced with x- I n taking the 
N — > 00 limit there are terms that diverge linearly with N and terms that behave as log x 2 . 



1(3 



Respectively these are 

[x r{vi— +v{ h Nvfw! — J 

m X (55) 

loga; 2 (x^i(l + J^) + + 3Nv 2 Wl X 2 ) . 

Setting these to zero fixes mi, Xi as i n with A replaced by x- We choose to keep 

the A shift, as is common in the literature [El E3 El 120] ■ If one wishes to translate our 
results, in the N — > oo formalism, only a simple substitution of A — > Xi/N is needed in 
the values of the diagrams given in Appendix O What is left is the finite correction to the 
leading Dyson equation for the \P propagator: 

r( A s ) + u + v 2 ^[ + Nwv 2 & 2 = . (56) 

Following the same procedure for the S, a, and u Dyson equations, one gets the finite 
equations 

p( A s ) +w + v + w 2 E[ - v 2 Yl 2 + Nw 3 E' 3 - 2Nv 2 wZ' i = 

P (A a ) + ^w + jw 2 n[ + ^w 3 W 2 - ^-wv 2 W 3 = (57) 
r(A ff - 1) + Nv + Nv 2 F[ + N 2 wv 2 F 2 = , 

where E-, IT-, and F[ are the finite parts of Ej, IT, and listed in Appendix O The 
minus signs in ()57|) come from fermion loops. From each Dyson equation one also gets a new 
determination of mi and Xi- 

1 w 2 S 1 - v 2 S 2 + w 3 S 3 - 2w 1 v 2 S i -w 2 S 1 + v 2 S 2 - 2w 3 S 3 + A Wl v 2 S 4 

m i = ic ; Xi ~ 



A wi + v\ W\ + V\ 

mi = -^K^i + w\P 2 - vjP 3 ) xi = -wiPx - 2wjP 2 + 2vjP 3 (58) 
1 . 

mi = - — (fiYi + WxViY 2 ) Xi = ~v\Yi - 2wxv x Y 2 , 

where P iy Si, Yi are the residues of Ilj, Ej, and F iy respectively 

Fortunately, the four seemingly independent determinations of mi and xi & U agree, as one 
can check by explicitly evaluating ()52|). (|54|1. and (fo"S|) using expressions from Appendix O 
with As = fi — 1 and A CT = 2. This provides a check that the renormalization procedure 
we have chosen to cancel the divergences of higher-loop graphs is consistent. Other schemes 
change the values for individual amplitudes, but the critical exponents remain the same [To] . 
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Interestingly, there is yet another consistency check on \i- One can show from (|54|) or 
dHHD that 

Xi = /i(2/i - 3)77a . (59) 



This is seen to comply with a scaling law formulated for the bosonic model in 

2A = 2/i - A a . (60) 

That this relation is also valid in our case can be seen by applying the Callan-Symanzik 
equation near the critical point for (a(p)a(—p)) or (u(p)u(-p)) propagator. 
Now we can solve (|56 )) -([57 |) by eliminating w and v: 

r(A*) = ^r(A CT - 1) + ±{F[ - + ^t(F> - d>' 2 ) . (61) 
Expanding Ag, A CT in l)fil|) . we can determine 772 



r/2 J- 
77? " 2/i 
/i 1 / 



f| = — + (/i - 1)(2/X - 1) (-1 + 7T COt jMK + 5(2/1 - 2)) 



/i-2 2(/i-T 



1 - ^ - 2) (^5(2/, - 3) - 5(/i - 1) - j±j + — - 2) 



(62) 



where 5(:r) = i()(x + 1) — ip(l) and the -B(x) function is defined in the appendix. The first 
line just comes from the Hatree-Fock diagrams, i.e. by iterating the 1/N Dyson equation to 
the next order, the second is the contribution of $1, Fx and the third comes from $2, F%. We 
also can determine the values of u>2, v 2 as 

^ (3 - n + {fi - 2)vr cot /itt + (/i — 2)5(2// - 3)) 



- fi ((7/i - 9)B(jjl - 1) + (13 - 10/1)5(2) + (3/i - 4)5(2// - 3)) (63) 

fi(fi- 1) 



fj, + 2/i(/i - r 



2/i -3 



-Jl- = 3|-_-£ + /,(/, - 2) + 2/i(/i - 2) (5(2) - B(n - 1)) . (64) 

TfrVx T]( 2/i 2 

In (|63|). the first three terms come from iteration of the first order equations, while in (|64[). 
the first two terms come from such iteration. 
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4.2 Calculation of A2 

As in section 13.21 the calculation of A2 through order 1 / D 2 requires evaluating each graph 
with one propagator altered from its leading power behavior (e.g. Tss/x 2As for an S propa- 
gator) to its sub-leading power behavior (e.g. r ss^'ss/ x 2As ~ 2X ) ■ The four Dyson equations 
lead to four linear equations in the quantities T' ss , T'^, T' aa , and T' uu : 



where for example we denote by all the diagrams that appear in the S propagator where 
the \l/ propagator is corrected. As in section E~T1 the amplitudes diverge when 2A5 + A a — 
2/i — > 0, and the same procedure described there to regulate and subtract the divergences 
and to remove terms proportional to log x 2 carries over to the present case. The finite parts 
of all the quantities in (|fj5|) are given in Appendix [0 as well as some further remarks on 
their evaluation. 

The system f!65|) must have a nonzero solution for the T's, so the determinant must be 
zero. This determines A2. A way to calculate the determinant to sufficient accuracy is to note 
that r(As)s(As, A) ~ 1/^V , and then expand the determinant into three 3x3 determinants, 
i.e. expanding in the line of \1/ field Dyson equation. All terms have to be expanded up to 
1/N 2 accuracy. One also notes that A2 only appears in the expansion of p(As)q(As, A) at 
this order. So A2 is going to be a linear combination of the various sums of diagrams given in 
the appendix, factors of u>2 and t>2, and terms that come from iterating the 1/N equations. 












0. 







(65) 
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The final result is quite involved and we prefer to give it implicitly as 

X > 1 25 ■25-gg±4- , 5 r -jtZj + 50Q.-l) 



r/ 2 2(^-1)2 2(/t-l) 2(/i-2) 4(2// -3) 2 

, / 19, N 45, n9 3, A /i 2 (2/i-3) 2 
+ (A* - 1 —=-ilM - 2) + -(/i - 2) 2 --(2n-3y 



2 Kn ' 2 V ^ 7 8(/i- 1 
2) "75(2/, -3) 



4 4^ y (A* -2) (2/1-3), 



- 2/x(/t - l)€f 2 - (| + 2/x(2/x - 3)) * 2 - 2/t(/x - l)(/t - 2): 
-M.-l) (2 ^ 3)2 ^-2^1^-/,^-2^ S < 



(66) 



771^1 



At -2 ?7iwi (/i- 2) 2 7/1V1 V/i-2 / r/i?Ji 
F 3 

-3/t(2/i-3)(/i- 1)— (vTcot/i7r + #(2/i-4)) . 

TJiVi 2 

The right hand side is a function of \x which can be obtained explicitly by substituting the 
expressions (J33J), @3J, (HJ), (|I27jl . {HID, (EH), and ffT35|) into fU). 



4.3 Calculation of the beta function 

As explained in section I3~T1 we can calculate the beta function once we know A. Noting that 
(166)1 gives the 1/iV 2 expansion term and subtracting 2Ai we find 

A = e/2, AxOO = ^ + £ - ^ + C(e 5 ) (67) 



A 2 = + |C(3)) 6 4 + 0(6 5 ) . (68) 

Using (|29j). we compute the beta function for the heterotic string in a constant curvature 
background: 

l„j g'D (, . D\ 9 5 B ! (i D\ 3„„, 5n2 , „, 1 



= -zv - ^ - v y i + 2 j - v u - jJ - 2 c(3 >* D + ■ (69) 

It is obvious that there is agreement with the first two loops of the expression where we 



use 



fl(o) = m^- = —g ij 0a (70) 

and (J2J). We do not know of any calculation of the beta function of the heterotic string in 
the minimal subtraction scheme that goes beyond two loops. In [2H 122] the beta function 
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was computed using the background field method, and found to be in three loops 

^ f—'R pkmnp pi \_ p p Inp p m ^p p/mnp p A (71^ 

Hij 8 \ 2 mnp ^ l m i 3 n P ^ jl ^imnp J ■ \ / 

The appearance of the Ricci tensor means that it is not minimal subtraction. Divergences 
involving the Ricci tensor can only appear through closed loops where at least one propagator 
starts and ends at the same vertex. Within the minimal subtraction scheme, at more than 
one loop these terms combined with their counterterms never produce a simple pole [2*3| I24j . 
A small check of our result comes from the famous C(3) term, ^^R mhki R jrt m (R k qs T R tqsh + 
R k qs t R hrsq ). This term is identical in the bosonic |231 I2H1 EI]- supersymmetric [231 121]; 
and heterotic [2H] cases. In an expansion of the Virasoro amplitude, it is associated with 
the constant term in an expansion in the Mandelstam variables s, t, and u. At loop order 
n + 1 in NLcrM calculations, it seems likely that the coefficient of ((n) is the same for the 
bosonic, supersymmetric, and heterotic cases (see fH] for a comparison of the bosonic and 
supersymmetric cases). 

In |29( |2H ED]; the absence of a three-loop term of the form a' 3 R 3 was noted. The three- 
point scattering amplitudes suggest that there are also no RF 2 or F 3 in the effective action. 
One knows that identifying the gauge connection with the spin connection in the heterotic 
string effective action will give the superstring effective action, where there is no a' 3 term. 
So if there were any R 3 terms in the heterotic case it would not be possible to cancel them. 
All this seems in conflict with the (|69[). where the term proportional to g 4 would seem to 
correspond to an R 3 term in the effective action. But it should be noted that the relation 
between the effective action and the beta function is [3^ 

2 «V|^ = Kp kl , (72) 

where KfJ can be computed perturbatively. In the bosonic case, this was done in the minimal 
subtraction scheme in [3D 132]; i n the heterotic case, this was done in a different scheme in 
|211 122] ; but we do not know of a minimal subtraction calculation of K\ l - in the heterotic 
case. In the bosonic case, K k J receives contributions starting at two loops, and it can be 
shown that this is compatible with an independent calculation of the effective action using 
scattering amplitudes. The same thing may happen in the heterotic case: in particular, R 3 
terms could indeed be absent from S e g, and the g A term in (J69|) could come entirely from 
K k J. A similar conclusion is reached in |3"3"1 IM] where it is shown that the beta function of 
the heterotic string in the presence of background gauge fields has a term at three loops that 
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behaves as F 3 , even though no corresponding term is present in the the effective action. 

Finally, it is possible to make a statement about the three-loop structure of the beta 
function in the a' expansion. Excluding the Ricci tensor and the Ricci scalar, since the beta 
function is computed within the minimal subtraction scheme, the terms that are third order 
in the Riemann tensor and are compatible with the g A terms in are given by 

I D D rnlr p k ^ p pfcmnr pi \ (7Q\ 

U J ' lr ~ ~\Q ik[ i ft mnr J ■ \(<J) 

4.4 Singularities of the critical exponents; central charge of the 
CFT 

Because A involves products of V functions it is natural to investigate the location of its 
singularities closest to the origin, as in [I]. Because A^ 6 * = ^X\ os , the location of the pole 
of Ai coincides with the pole in the bosonic case, with half the residue. One also has to 
note that behaves as r]i ~ — ^jjt^i * ,e- ^ nas a srm pl e P°l e at e = — 1. But Ai's first 

singularity is at e = —3, since the pole of rji is canceled by a similar pole of Xi- Examining 
term by term the structure of A2 it is easy to see that the singularities of A2 come from the 
r\\ factor that multiplies the whole expression and from the three-loop diagrams that 
have the lagrange multiplier field propagator corrected, i.e. n 2o -, IT^, F 2u , and F 2a - Since 

RM ~ j^- r ~ (74) 

and A2 has terms that behave as ~ R3V1 and R2Tjf times a fi polynomial with no zero at 
fi = 1/2, we see that it has a fourth order pole. In all, one finds 

Xl = Jjm + 0(1) 6l= _!^|+M +0(1) (75) 

^ = T^rV + °^ + i)" 3 ) h - = + ° {{K+ 1)_2) • (76) 

The l/(« + l) 2 term in 62 comes only from the factors R31H and from the Hartree-Fock 
diagrams. The singularities in the heterotic case are at the same locations and of the same 
order as in the bosonic case. 

Because of the sign of 62, there is clearly a zero of (3(g) (computed through order 1/-D 2 ) 
for negative g, close to k — — 1. The same caveats discussed in P apply: higher order terms 
in 1/D could conceivably cause this zero to disappear or move significantly. In section [HI we 
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will comment further on higher-order corrections. For the remainder of this section we will 
assume that the computation of (3(g) that we have carried out is precise enough to describe 
the zero correctly. 

The zero of j3(g) arises through competition between the one- loop term (corresponding 
to Einstein gravity) and 62 (corresponding to a combination of all a' corrections to Einstein 
gravity). Because the geometry has string scale curvatures (more precisely, L 2 ~ Da') 
there is no reason to think that the worldsheet central charges are particularly close to the 
flat-space results. Fortunately, one can calculate the central charges using Zamolodchikov's 
c-theorem: 

The result (|77|) holds for both the holomorphic and the anti-holomorphic sides: c and c 
differ by a constant. To derive the prefactor on the right hand side of 1)77)1 . one can consider 
two-point functions of the graviton perturbation Oy = ^^OXidXj + ^-^d^f around flat 
space, as is done in £Q. 3 This prefactor receives higher loop corrections, and knowing K^J 
in higher loops, one can in principal compute them. As in the bosonic and supersymmetric 
cases, the results suggest that with increasing D the critical point moves closer to k — — 1: 
integrating (|77)l leads to 



^ 3(£> + l) f Kc , I f 6i(k) b 2 (K)\ 3,„ 1W1 



(79) 



where we have noted that the central charge of the holomorphic side in flat space is c = D+l, 
while for the anti-holomorphic side it is c = + 1). The approximate equalities arise from 
dropping the &i(k) and 62 («0 terms from the integrand: their only role at this level of 
approximation is to set k c . As k c gets closer to —1 (i.e. as D becomes large), the central 
charges converge to 

c=|(D + l) c = 3(^ + 1). (80) 



3 Another way to derive (|77|l is to use the relation of the central charge to the spacetime effective action. 
At least up to two-loop order, the effective action at the fixed point is equal to —c/2n 2 D a' |35l 1361 1?7] . Also 
up to two loops, the K^J of Ij72(l is simply given by a product of Kronecker S's, as in the bosonic case |H2j . 
Using the fact that in symmetric spaces 



mm = -gPij 9 ij -^j=g§- g , (78) 



one indeed ends up with 1)77(1. 
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The result ()8U|) for c is the same as in the bosonic case, while for c it is the same as the 
type II case [TJ. As in jTj, (|5U|) appears to set only an approximate upper bound on the 
central charges. The dominant error in the calculation [7§] is from the uncertainty in the 
prefactor in (|77jl . Analogous to the speculations in £Q , it is conceivable that the expressions 
()80[) might in fact be exact. But this would require a significant conspiracy between the 
prefactor in (J77|) and the beta function. 

The fact that the location of the critical point at finite D is so close to the singularity of A 
means that the critical exponent A evaluated at the critical point is large and positive. This 
leads to an operator with a large and negative dimension, which appears to violate unitarity. 
However, one could hope that a consistent GSO projection would project this operator out 
of the spectrum. 

5 Discussion 

The existence of the AdSo+i critical point depends on competition between one-loop and 
1/D 2 effects. It would therefore be instructive to compute the beta function through order 
1/D 3 and see whether the fixed point persists. Given that the number of diagrams needed 
for the computation at the next order grows significantly, the shortest path seems to be 
calculating % 3 and using (|5U|) to deduce A3. However, note that for the calculation of xs 
one needs to derive the residues of diagrams at order 1/D 4 , which include some six-loop 
diagrams. 

There is some reason to think that the singularities of A at order 1/D 3 are no worse than 
at order 1/D 2 : examining the diagrams needed for the Dyson equations of the Lagrange 
multiplier fields, we see that at order 1 / D 3 these come from either inserting a o or u prop- 
agator in the 1/D 2 diagrams or inserting a loop of S or \1/ in the middle of the diagram. 
The computation for the diagrams that come from inserting a a or u propagator can easily 
be seen to be reduced to the sum of diagrams similar to TI2 or n3 with one different expo- 
nent. A naive calculation does not produce any worse singularities than the ones already 
contained in U 2 and n 3 . However, one also has to compute the more difficult diagrams with 
the additional S or \1/ loop. 

It is evident that the methods of [E1E3> has many advantages over calculating Feynman 
diagrams in momentum space. In the latter approach one encounters difficulties already in 
calculating second order diagrams, since the propagators of the Lagrange multiplier fields 
are in general hypergeometric functions. It is noteworthy that even though we start from 
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d = 2+e dimensions, one can calculate critical exponents of the O(N) model in any dimension 
2 < d < 4, and there is agreement with the results in three dimensions jSH] i n the bosonic 
case. 

Perhaps the methods of [HI Ei] could be applied to a related quantum field theory: 

S = j d d x Qv$ ■ V$ + ■ $ - , (81) 

which for d — 3 is the proposed dual of an AdS^ vacuum of a theory with arbitrarily high 
spin gauge fields [SH] • What makes (jHTj) susceptible to a position-space treatment analogous 
to those in jHllIS] is that only cubic vertices are involved. It would be interesting to compute, 
for example, the four point function of a to order 1/N 2 and compare it to the corresponding 
Ad S4 calculation, as is done for example in [ID] at order 1/N. 
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A Anomalies 

Since we have coupled only the right moving fermions to gravity it is natural to investigate 
whether there are anomalies. These are related to a breakdown of general coordinate invari- 
ance or local Lorentz invariance. We will investigate only the latter, as is usually done |41j . 
Indeed when there is a coordinate anomaly one can add a counterterm to the action and 
convert it to a Lorentz anomaly [12] • It is convenient to use the tetrad formalism. Local 
Lorentz transformations in this formalism are given by 

e'/(x) = e^(x)Q%(x) . (82) 

The Riemann tensor can be written R^q, with mixed spacetime and tangent space indices, 
and can be regarded as a two form R 2 . We only have to worry about the massless fields 
of the supergravity sector [3J. The anomaly polynomials for the spinor and the gravitino 
contain only terms that are proportional to polynomials in Tri?2 2m - In the AdS space that 
we are interested in we can calculate 

Rfj,v a bRn\ C a = Y2(5 C kR[iiy\b + & 'xR^ubn) (83) 
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a + ft — [i 



a + — fi 



a + j3 — \l — 1 



u(a, f3, 2/x — a — j3) 



(o-l)(2/i-a-j9) 



v{a-l,f3-l,2[j,-a-l3+l) 
2(a-l)G8-l) 



Figure 6: Products of two propagators are related to a single propagator. So, in a two-loop 
diagram, (Ex for example) inserting a point in one of the three propagators that connect to 
one internal vertex can make this vertex unique. The dotted lines denote fermions. Since 
we are dealing with chiral fermions, taking the trace in the third graph only produces one 
half the full result, u is equal to v(xi, x 2 ,x 3 ) = it Yl^=i a ( x i)- 





7 v(a,(3,-y) 



H — a 

Figure 7: An identity that allows the integration of a unique vertex. Only bosonic propaga- 
tors are shown. Similar identities with fermions can which are used in our calculation can 
be found in 



which when antisymmetrizing to get the wedge product returns zero. So R 2 A R 2 = in our 
case, and we do not have to worry about the gravitational anomalies. Another way to view 
this is to say that the field strength H3 = dB 2 obeys the modified Bianchi identity 



dm 



1 

Air 



(Tr R 2 A R 2 — Tr F 2 A F 2 ) . 



34) 



A three form H3 obeying (jbljl is required for cancelation of perturbative heterotic string 
worldsheet anomalies, as is briefly reviewed in [32] ■ Because Tr R 2 A R 2 = 0, this is trivially 
satisfied. 
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B Position space methods for calculating graphs 



There are 11 diagrams needed for the calculation of t] at 1/D 2 . We designate them by S, 
II, <£>, F. The way to compute them was developed in EI] for the bosonic graphs and 
extended to include fermionic graphs in The main advantage of the method is that 

there is no need to explicitly evaluate any Feynman diagram. Here we will only give a few key 
observations that facilitate the evaluation of the diagrams. The first observation is that the 
chain of two propagators is equal to a propagator times a prefactor. Graphically this is shown 
in Fig El where v(xi,X2,X3) = ir^ Yli=\ a i x i)- The third exponent x 3 is determined by the 
"uniqueness" ^3] requirement Y2i x t = for the bosonic graphs and x« = 2/i — 1 if there 
are one or more fermion lines in the graph. An identity exists for a three point vertex, which 
is similarly related to a "unique" triangle, where now the uniqueness requirement is that 
Y2i x i = f 1 - ^ there are one or more fermion lines the uniqueness changes to ^2 { x i = f 1 + •"•> 
and the results of jTB] are unchanged in our case. There is no similar identity for a four 
point function, and for the (1,1) supersymmetric model that means that one has to retain 
the auxiliary field F. In computing the values at order \ is set to zero [T5] . For a non- 

zero A the diagrams, for example the self energy of \& designated A, lose their uniqueness. 
However one can subtract from A a graph B that has the same divergent substructure as 
A, but can be calculated for an arbitrary A. So one has to compute (A — B) + B. Since B 
can be calculated for arbitrary A and contains the divergence, one can evaluate (A — B) at 
zero A when both diagrams become unique. A valuable first step in the calculation is the 
evaluation of all the self-energy graphs, which we will not include here since it was done in 
detail in ^SJEl- We just note that the most basic tool is the insertion of a point facilitated 
by the fact that we can write a propagator as a product of two. Then one can choose one 
of the exponents in such a way that the vertex that the propagator is attached to, becomes 
unique. 



C Calculation of the graphs needed for 772 

We give the results for the various graphs occurring in the 1/D 2 calculation. We only give 
the simple pole term and the constant term in an expansion in A. The purely bosonic graphs 
were calculated in [T^]. Compared to the calculation of the fermionic graphs in [20], there are 
differences that have to do with taking the trace of fermion loops, i.e. some factors of two in 
bosonic diagrams with fermion loops. Otherwise the calculation is almost identical. We find 
small discrepancies with [20] in some of the diagrams, mostly factors of 2 and some minus 
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signs. The most notable difference is $1, where the residue has a different denominator. We 
believe that our value is correct, since it leads to the same \i as the evaluations from the 
other equations ([5%]). 

B(x) =ip(x)+tp(fi-x) (85) 



2 7 r 2 ^ 2 (A 5 )q(A tT — 1) / A 
Aa(A a - l)r(/i) V 2 



B{A a -l)-B{A s ) + 



A„ - 1 A, 



Sa = ^(A s)a3 (A : H, + A s -A,) Q + A [B(K) B(As)] 
^ 7T 4 ^ 3 (A S )« 2 (A (7 - l)a(A CT )a(A s + // - A a ) 



AA 5 (A 5 + / i-A ff )(A ff -l)2r(/x) 



A 

- + 2 



B{A a ) + 35(A CT - 1) - 4£(A 5 ) 



A„ - 1 A, 



(86) 
(87) 



(89) 



(90) 



n 2 



vr 4 ^ 3 (A s )a 3 (A CT )«(A 5 + // - A a ) 



n. 



(1 + A [4B(A ff ) - 35(A S ) - 5(/i + A 5 - A,)]) 

(91) 



Ar(Ai) 

7r 4 ^ 3 (A s )a 2 (A CT - l)q(A s + //-Aj 
2A«(/i - A ff )A s (A s + /i - Aq-) (A(j - l) 2 r(/i) 

x (\ + A [2B(A„) - 3fi(A 5 ) + 2i?(A cr - 1) - £?(A 5 + ^ - A CT ) 

1 2 1 

+ 



A s A CT 1 A 5 - A CT + /i 
vr 2 ^ 2 (A 5 - l)a(A ff 



AA S (A S - 1)I» 



1 + A 



B(A a ) — B(A$ — 1) 



x 1 + A 



7r 4 ^ 3 (A s )a 2 (A CT - l)q(A g )a(Ag + ^-Aj 
AA 5 (A 5 + /U -A a )(A (J -l) 2 r(^) 

5(A CT ) - 25(A 5 ) + B(A a - 1) 



A, 



(92) 



(93) 



(94) 
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Fx 



2^a 2 {A s )a{A a - i; 
AA 5 (A (T - l)T(n) 



1 + A 



B{A a - 1) - B(A< 



2A< 



+ 



A CT 



(95) 



x 1 + A 



1 

A7 



7r 4 ^ 3 (A 5 )g(A CT )a 2 (A fT - l)g(A s + /i - A„) 
AA 5 (A s + /i-A (T )(A (T -l) 2 r( / u) 

£(A CT ) + 3£?(A CT - 1) - £?(A 5 + ^ - A CT ) - 3£(A S ) 

3 1 



(96) 



A CT -1 A 5 + /i-A a _ 
Note that there is a similar three- loop diagram with n 3 with the role of the u propagators 
interchanged. But it scales as 1/N 3 . A very useful identity for the evaluation of various 
quantities given in the text is B(x) = B(fi — x). 



D Calculation of the graphs need for A2 

In this section we give the formal expressions for the sums of the corrected diagrams 
(|1(J7|) . the values of the 43 individual diagrams that contribute (jl08|) - (|126|) . and finally the 
explicit form of the sums ()127j) - (jl35p . Beforehand one has to define the functions appearing 

as 



R x = - 1) - 

R 2 = t//(2fji - 3) - f{2 -fj)- $'{n - 1) + ^'(1) 
R 3 = tP(2f, - 3) + V(2 -y)- ^ - 1) - ^(1) , 

where if}(z) = dlogT(z)/dz. 

E s = 2w 2 Y, 1Sa + w 2 Z 1Sb - v 2 Z 2S +Nw 3 (2i; 3Sa + 2E 35b + S 35c ) 

-2Nv 2 w(T, 4S a + S 4 56 + S 4Sc ) 



(97) 



(91 



= — 2v 2 T,2^ — 2Nv 2 w(E^ a + £4^) 



(99) 



2w 2 E la + iW(2£ 3ffa + £ 3(j6 ) - 2Nv 2 wY> 



lcr 



(100) 



-2v 2 E 2u - 2Nv 2 w(T, 4ua + E 4ub ) 



(101) 
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(102) 



$ 5 = 2vH ls + 2Nv 2 w(^ 1Sa + $256) $ u = 2t; 2 $ lu + A^ 2 w$ 2u (103) 

n 5 = 4u; 2 n 15 + Nw 3 (AU 2Sa + 2U 2Sb ) - 4i\Wn 35 (104) 

n CT = u; 2 ni ff + 2Nw 3 U 2r7 = -2i\Wn 3 * U u = -2Nwv 2 U 3u (105) 

i 7 * = 2v 2 F^ + 2iVt; 2 W 2F 2 ^ F 5 = 2v 2 F ls + 2iYVu>(F 25a + F 2Sb ) (106) 



F u = w 2 F lu + Nv 2 wF 2u F a = Nv 2 wF 2a (107) 

There are 19 diagrams associated with the S-propagator and eight for the propagator of the 
other fields. The value of each graph of each graph is given below for completeness. The 
bosonic ones were calculated in ^3], while similar to the fermionic ones were done in |19j . 



7T 2 ^ „ 7T 2 ^ 



nSa — 7 n , r9 , x ^>lSb 



(/x-2)TV) ° (A* - 2) a F»(/i) 

2(/j+l)^ _ (/i 2 - 3/1+1)1(1-/1)^ 

J2S " M/i-i)r 2 (/i) 35a ( / i-2)3r( / i)r(2 / i-3) 

7r 4 T(2-/i) / D 2^-3 



108) 



(109) 



(2-//)r(//-i)r(2//-2) V (/x-2) 



3i? i + 77^ o\2 ) ( 11Q ) 



^r(4-^) gTr^rq-AQ , . 

35c (/i - 2) 3 r( Ai - i)r(2/i - 4) 4Sa r(/x)r(2/x - 2) 1 j 



7r 4 T(2-/x) / 2/x-3 



r(^)r(2^-2)V (/i -i)(// -2 



3i2i + t ^ (112) 



_ 2(/i-3)(2/i-3)r(l-/i)7r^ _ - 

" (2-^(2^-2)^) Ew " " (U3) 
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u * a ~ 2(/i-2)r(/i)r(2/i-2) " ~ $2Sa L4 * fe " 2(2-/i)r(/i)r(2/i-2) ~ ~ $2S6 

(114) 



tt^V-Sm+I) 7r^(2 ; u 2 -7 / i + 4)r(l-^) 



, 37T^r(3 - /X)^ _ 7T^(2/X-5)r(l-/x) -j-j- , . 

J3CTfe " (2-/x)3r(/i-i)r(2^-2) " 11256 L4CT " (^-2)r(/x)r(2^-2) " 1135 (116) 

y y ^(4/x-9)r(l-/.) 

2u r 2 (/i) 4ua 2(/i-2)r(/i)r(2/i-2) " 2sa 1 J 



37r^r(2- / i)/? 1 

J4ufe 2(/x - 2)r(/i)r(2/i - 2) 256 



(118) 



$1 "=(^3%) (U9) 

7r^r(l-^)(V-llm + 5) gTT^gi 
2u 2(1 -/,)(/, -2)r( A *)r(2 A * -2) lCT (2 -/*)(2/*- 3)1*^-1) 1 J 



2(2 - /i) 3 r(/i - l)I\2/i - 1) V (/x- l)(3/u- 2) v n-2- 



(121) 



n r(2-/i)7r^ / 2 2i? 3 (^-2)-2 \ _ 

^ " "4(,x- 2)^)1(2/, -2) l^ 1 " R2 ~ R * + ( Ai -l)(2 M -3) J " (122) 

" (2^-3)1(^-1) " $2CT " (l-ffl (123) 

7r^(4/i 2 -ll^ + 5)r(l-^) 

2 * 2(l- / ,)(/x-2)r( / i)r(2/i-2) 1 ] 
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7T 



2/< 



'15 



2u 



r 2 M 

7r 4 T(2 - fj) 



Flu 



40i-i)V-2) 2 r(Ai)r(2/i-2) 

, 2(^-2)^3 



3(/i - 1)tt 2 ^ 



(A« — l)(6J?i — ife - 



+ 



4 



2// — 3 2/x — 3 // — 1 



(125) 



(126) 



D.l Summing up graphs 



In this subsection we give the explicit values for the various sums appearing in the A2 
calculation. We have omitted the 1/N 2 factor that multiplies all of the diagrams so that the 
expression (jo^j) for A2, does not contain any factors of N. 



•1 + H±| 4 



- + 2/i(2/i - 5) + 2fi(2fi - 3)(/x - 3)(2 - (ji - 2) 2 ) 



)JL — 1 (/! — 1 

- 4/i(/i - 2) - 6/i(// - 2)i?x 



/i 2 (/i-2)(2/i-3) | 
/i — 1 

/i(3 + /i(3/i - 7)) 



Z(ji-l)(ji-Z)R x ) =$5 



(127) 

(128) 
(129) 



= r/iui(/x(8 - 4/i + 3(// - l)i*i) = 



(130) 



3//(/4 -!)(//- 2)^i 
^ = ^1 — - = (n*)/2 



2/i-3 



$u -1 + 6^-4/x 2 



771^1 



n CT f3fi(ji-2) 



r/ivi V 2 (2/^-3) 



Hi 



At-2 



6i?i — i?2 — Ho 



6i2i -R 2 -R 



rjivi 



2 2(M^-1) 
(2/1-3)^-1) 



, (//-2)iZ 3 -l 
'(2/i-3)(/i-l) 



2F CT 



(131) 
(132) 
(133) 
(134) 
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3 

4~ 



jj-2 



i 



6R 1 



i? 2 + R3 



2 2(^-2; 
(/i-1) 2 ' (2/x-3)(/i- 1) + 2/i-3 



+ 



(135) 
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